A chiral random matrix model with 2+1 flavors at finite temperature and density * 
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Phase diagram of a chiral random matrix model with the degenerate ud quarks and the s quark 
at finite temperature and density is presented. The model exhibits a first-order transition at finite 
temperature for three massless fiavors, owing to the Uyi(l) breaking determinant term. We study 
the order of the transition with changing the quark masses and the quark chemical potential, and 
show that the first-order transition region expands as the chemical potential increases. We also 
discuss the behavior of the meson masses and the susceptibilities near the critical point. 



Introduction. — Study of the QCD critical point (CP) 
P, [2, Q is an intriguing fundamental issue since its ex- 
perimental confirmation will yield a strong evidence for 
the QCD phase transition, and energy-scan experiments 
searching for the QCD-CP are being performed at Rel- 
ativistic Heavy Ion Collider at BNL. Although the exis- 
tence of the critical point in the QCD phase diagram is 
presumably inferred from the model studies and lattice 
QCD results0,i,[3, its absence is also a possibility[l,0]. 
In this paper, we adopt as a schematic model for QCD 
the chiral random matrix (ChRM) model which incor- 
porates the UA(l)-breaking determinant term0,[l]. We 
report the phase diagram of this model with the degener- 
ate ud-quark mass m^^ and the s-quark mass rus at finite 
temperature T and quark chemical potential ^. 

The ChRM models have been successfully applied for 
qualitative study of chiral properties of QCD[^, [l^. In 
a ChRM model the Dirac operator on gluon field back- 
ground is modeled by a matrix D in the space of constant 
modes with small Dirac eigenvalues, retaining the chiral 
symmetry {Z?,75} = 0. The partition function of the 
model is given as an average of det D over random en- 
semble of matrix elements, which mimics the complexity 
of the gluon dynamics. The finite T and /i effects arc 
treated schematically as non-random external parame- 
ters appearing in D. In Ref. [To| . the phase diagram of 
the ChRM model has been explored in the T-^ plane and 
a tri-critical point (TCP) is found on the phase bound- 
ary in the massless limit. The TCP changes to a simple 
CP when the quark mass is nonzero. This result is con- 
sistent with the phase structure obtained in other model 
studies with two quark flavors [l|, [3, [HI implying the sce- 
nario that the CP exists in the QCD phase diagram as 
an endpoint of the first-order phase boundary. 

Nature of the chiral transition in QCD is sensitive 
to the number of light quark flavors, especially to the 
value of the s-quark mass TOs- Unfortunately, however, 
the phase structure of the conventional ChRM modelfioj 
is independent of the number of flavors Nf. In order 
to remedy this problem, we have recently incorporated 



the UA(l)-breaking determinant interaction 0, in the 
ChRM model @ by extending the zero-mode space [12 
with the instanton gas model picture in mind. This is 
the first ChRM model which describes the Nf depen- 
dence of the chiral transition allowing us to explore the 
phase structure varying the parameters mud and rUs in 
addition to T and ^. 

Model with determinant interaction. — The chiral sym- 
metry breaking manifests itself in the nonzero density of 
the zero Dirac eigenvalues through the Banks-Casher re- 
lation [l^. The origin of small Dirac eigenvalues may 
be instanton configurations of background gauge field 
and other nonperturbative gluon dynamics. Here in our 
model we divide these fermionic modes into two cate- 
gories, N+ and iV„ topological zero modes associated 
with A^-i- instantons and A^- anti-instantons respectively, 
and 2N near-zero modes generated by other complex 
dynamics @, |T2|- The Dirac operator D then approx- 
imated with a martix of 2N + + iV_ dimensions 
with Af-|_, 7V_ and N being of the order of the space- 
time volume 0{V). The thermodynamic limit is taken 
a.s 2N + N+ + —> oo. Note that N± should vary 
depending on the instanton distribution. 

For fixed number of zero modes the model partition 
function is written in the chiral basis as 

dR e-^s'*^««' J] detiD + m/), (1) 
/=i 



with 



D 



iR + C 
ii?t + C"^ 



(2) 



where R <E c''N+n+)x{n+N-) -g ^ random matrix fol- 
lowing a Gaussian ensemble distribution with the vari- 
ance l/(iVS]2) and C £ ciN+N+)x{N+N_) matrix 

representing the effects of T and /i. The matrix D has 
|iV+ — I exact zero eigenvalues when R and C are rect- 
angular, which is interpreted as a realization of the index 
theorem in the ChRM model. We adopt here the simplest 
form for C^: 
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where T and /i are schematic representation for the 
temperature and chemical potential effects, respectively. 
Note that D with 7^ is non-Hermitian whereas the 
partition function ([1]) is still invariant under fi —fj.. 
One should appreciate that the A^_|. x N- right-bottom 
block in C corresponding to the topological zero modes 
is set to zero. This seems a reasonable assumption if one 
notice that the finite T and fi effects are introduced as 
a boundary condition in the Matsubara formalism and 
that the localized topological zero modes will be insen- 
sitive to the boundary. This discrimination is important 
indeed in reproducing the physical T dependence of the 
topological susceptibility 

The complete partition function of the model is ob- 
tained by summing over and with a distribution 
function P{N±): 



7RM 



J2 P{N+)P{N^)Z, 



(4) 



N+,N- 



The P{N±) reflects a modeling of the instanton distribu- 
tion in the QCD ground state. The authors of Ref. [T3| 
adopted the Poisson distribution, which involves arbi- 
trarily large number for N± and results in a model with 
no stable ground state. Inspired by the lattice gas model 
within a finite box, we instead choose the binomial dis- 
tribtion 



P(iV±) 



p^±(l-p) 



lN-N± 



(5) 



where 7 is a parameter of 0{V'^) and p is interpreted as 
the probability for a single instanton to occupy a unit 
volume V/{'yN). This distribution sets an upper bound 
'jN for N± and gives rise to a stable effective potential 
as a function of order parameters. In fact, applying the 
standard bosonization procedure to (|3]), wc find 



^ dS< 



X det f [{S + M)iS^ + M^) - (a* + iTf] 
X det * [{S + M){S^ + M^) -{fi- ir)2] 
X [adct{S + M) + 1] ''^ [a det(5't +M^) + l]'''^ 
dS e-2^^^(^^ , (6) 



where we defined the effective potential i^{S; T, fi) in 
the last line. S £ C^f^^' is the order parameter matrix, 
and A4 is the mass matrix. The parameter a ~ p/{l ~ 
p). Note that the integrand of Z^'^ is a polynomial of 
S except for the exponential factor originating from the 
Gaussian ensemble distribution. Large values of S are 
suppressed by this Gaussian weight. 

The determinant term with the coefficient a represents 
the anomaly which breaks explicitly the U^(l) symmetry 
of the effective potential Vl{S) even when = 0. For 



S = (t)lN; (0 e K) with M = Q.,U simpfifies to 

(7) 



7ln|a0^^' -I- 1| 



We see that the anomaly term yields —aj(j)^f when 
expanded. In Ref. [1] we studied this ChRM model 
with two and three equal-mass flavors at finite T with 
/i = 0, to show a second- (first-) order phase transition 
for Nf = 2(3). 

Phase diagram and meson masses with 2+1 flavors. — 
Choosing S — diag(0iid, 0ud, 0s) in the 2-|-l flavor case 
with A4 = diag(mud7 Wud, ^s), we have 



- 7ln|aVud¥'s + 1| , 



iTf] 



(T- 
(T 



-T)) 
-T))] 



(8) 



where (ySud = (f'ud + "^ud and ips ~ (t^s + tt^s- The ground 
state is determined by the saddle-point condition 



an 

90ud 



an 

0, T7- 



0, 



(9) 



which becomes exact in the thermodynamic limit. 

Prior to the numerical analysis, we comment on the 
model parameters E,q; and 7. Setting E = 1 by redefi- 
nition of S and other parameters, we searched such a set 
of parameters a, 7, mud and mg that reproduces quanti- 
tatively the (ratios of the) meson masses in the vacuum, 
but it was unsuccessful. However, the model can describe 
the mass hierarchy qualitatively as seen below in Fig. [3l 
and we wish to study the model phase diagram as an 
schematic model for QCD. Note that all the quantities 
arc dimensionless in this work. 

We also remark that the anomaly term makes a 
symmetry-broken phase more stable. Indeed, no symme- 
try restoration occurs at finite T for aj > S^(= 1) with 
Nf = 2, and the situation is similar even for Nf = 3. 
Hence one must assume Q!7 < 1 for study of the chiral 
restoration. 

Finally, the transition at finite /i ^ with T = is 
first-order. It is seen in the simple case ([7]) because the 
symmetric phase (/> = and the broken phase (jj > fi 
are separated with the point (j) = where = 00 or the 
integrand of vanishes. This feature survives in more 
general cases with 2-1-1 flavors. 

Let us study the phase diagram in the T-m^^-ms space. 
Since our model shows a first-order transition at finite T 
for ?Tiu(i = "nis = and a crossover for large rriud and 
TOs@1j there must be a line of a second-order transition 
separating these two regions in the m-a^-ms plane. This 
critical line is determined by the condition fi^") = (n = 
1, 2, 3) with = where ^2 is a function of 

a single order parameter (f)-^^ with (ps eliminated by the 
second equation in 0. Note that fi'-^-' = is equivalent 
to det a^n{S) / a(j)ia(j)j = O {i,j = ud,s), which implies 
the vanishing a mass (see below). 
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FIG. 1: The critical curves on the mud-rris plane (a) for 
a = 0.3, 0.4, 0.5 with 7 = 1 and (b) for 7 = 0.9, 1.0, 1.1 with 
a = 0.5. The finite-T transition is first-order in the smaller- 
mass region and crossover in the larger-mass region. On the 
rris axis with rriud = 0, there is the Nf = 2 chiral symmetry. 
The TCP is denoted by a dot for each parameter. 




FIG. 2: Critical surface in the mud-ms-fi space with the 
parameters 7=1 and a = 0.5. A series of TCP's is denoted 
by a thick line, and the second-order transition with Nf = 2 
chiral symmetry occurs in the shaded area. 



We present in Fig. [T] the critical line projected onto 
the TOud-'T^s plane for several values of a and 7. When 
we increase the strength of the anomary term a and/or 7, 
the region of the first-order transition expands. For each 
parameter a TCP is found on the uis axis, where the 
Nf = 2 chiral transition changes from a first to a second- 
order one. Near the TCP the critical line behaves as 
(mj^^ — TTig) (X m^^^ as is expected from the mean-field 
Landau-Ginzburg analysis, which is clearly seen in Fig.[T] 
(b). On the other hand, the line smoothly intersects the 
mud axis with a finite slope^*! . In fact, the model with 
the anomaly term is symmetric under rriud ^ ^^^^-ud but 
asymmetric under rUs <-> —mg. 

Next we extend our calculation to the finite n case. 
Fig. [2] exhibits the phase diagram in the m^^-ms-iJ? 
space. We see that the region of the first-order transi- 
tion expands as /i is increased. This behavior indicates 
the existence of the CP in the T-^ plane with the phys- 
ical quark masses, provided that the finite-T transition 
at /U = is crossover. Varying the anomaly parameters 
a and 7, we have confirmed that the expansion of the 
first-order region with increasing is a robust result in 
our model as far as we keep a and 7 constant. 

The (screening) mass matrices for the scalar and 
pseudo-scalar mesons are defined as the curvature of the 
potential around S = diag((/)ud, <l>ud, i's)' 



M, 



ah 



S=S 



(10) 



s=s 



where 5" = ^a{o'a + i7i'a)/v 2 with real parameters da 
and tTq, and with Aa being the Gell-Mann matrices and 

Xq = y^diag(l, 1, 1). See Appendix for explicit expres- 
sions for M^i^. Because the quark mass term m^^ 7^ nis 



breaks the SU(3) flavor symmetry to cause nonzero mix- 
ing = 



channels, we diagonalize the matrices to get the mass 
eigenvalues corresponding to a and /o for the scalars and 
77 and 77' for the pseudo-scalars. One can show that the 
flat direction of il(0ud,0s) near the CP coinsides with 
the a fluctuation direction. 

In Fig. [HI we show the meson masses as a function of 
fi a,t T = Tc with parameters 7 = 1, a = 0.5, rriud = 
0.05 and mg = 1.0. The model reproduces the empirical 
hierarchy of the meson masses qualitatively in small fi 
region, thanks to the anomaly term. With increasing T 
and/or ^ the pseudo-scalar meson masses remain nearly 
constant while the scalar ones, especially the a mass, 
decrease. At the CP (T — Tcfi = Hc), the cr meson 
becomes massless. At higher fj. > fic, pairs of the masses 
M^-Mtt, Mk,-Mk, and Ms-M,f get almost degenerate, 
reflecting the approximate Nf = 2 chiral symmetry. 




^ = A/sn 7^ in both the scalar and pseudo-scalar 



FIG. 3: Mesonic masses as functions of /i with T = Tc fixed 
for parameters 7 = 1, a = 0.5, mud ~ 0.05 and rUs — 1.0. 
The critical point locates at {Tc,fic) = (0.817,0.738). Thin 
(thick) lines indicate (pseudo-) scalar mesons. 



We thank T. Hatsuda's comment on this point. 



Generalizing the mass matrix to = Xa{sa + ipa)/ 
with real sources Sa and pa, the scalar susceptibilities are 
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defined as 



d'^n{S{M);M) 
dsadsb 



(11) 



where S{A4) solves © for a fixed Ai, and the derivatives 
are evaluated at M ~ diag(mud, "iud, "ig). The pseudo- 
scalar susceptibilities x^l ^re defined similarly as the re- 
sponse to Pa- Wc find the (pscudo-)scalar susceptibilities 
Xab Xab ^ diagonal form 

xT = 6a,J:'{J:VM:r^'~i) (12) 

for a,b ~ 1, . . . , 7. There is an mixing term for a, 6 = 0, 8 
due to the SU(3) breaking: 



s(ps) 



Xab 



(13) 



which becomes diagonal when the mass matrix M'^^p*'^ 
is transformed to be diagonal. Note that the scalar sus- 
ceptibility in the a channel diverges when the (screening) 
mass M„ vanishes at the CP. So do the quark number sus- 
ceptibility Xq = —d'^^/dy? as well as the 'specific heat' 
XT = —d^Vl/dT^, through the mode mixing generated 
by the finite condensate cpud and ips ■ 

Conclusion. — As a schematic model for QCD, we have 
analyzed for the first time in the ChRM model the phase 
structure with 2-1-1 flavors at finite T and /x, which be- 
comes possible by the inclusion of the Ua(1) breaking 
termQ. We have drawn the critical curve separating the 
first-order transition region and the crossover region in 
the mud-'m-s plane with fi = 0, and we have shown that 
the first-order transition region expands as the strength 
of the anomaly term is increased. 

Extending the model to the finite fj, case, we have 
shown the critical surface in the TOud-'7is-/i space. We 
have found that the first-order transition region expands 
with increasing /i, which is a supportive result for the ex- 
istence of QCD-CP on the T-fj, plane: one encounters a 



CP as fj, is increased from zero, provided that the finite-T 
transition is crossover at /i = 0. The meson mass hierar- 
chy in the vacuum is qualitatively reproduced with the 
Ua(1) anomaly term and the SU(3) flavor breaking. At 
the CP the (screening) mass of a vanishes. Although we 
treat the model parameters independent of T and /i, pos- 
sible rapid quenching of these parameters at flnitc fi can 
give rise to a shrinkage of the first-order transition region 
jig . Furthermore, the role of the vector interaction be- 
tween the quarks [l3| deserves further study. These are 
beyond the scope of this schematic model. 

Appendix. — The mass matrix A/^^ is diagnal except 
for a,b = and 8. For a = 1, 2, 3, 



(14) 



and for a = 4, . . . , 7, 



M^^^ _S ± Re(^2^_,2)(^2_,2) ± ^^^^-PT ' 

(15) 

where ip^d = </>ud+TOud, '^s = (f's+rris and z = fj,+iT. The 
upper (lower) sign corresponds to the (pseudo-) scalar 
meson. The elements M^j, for a,b = 0, 8 are written 
concisely in another basis of X^d = diag(l, 1, 0) and Ag = 
diag(0, 0, V2), instead of Ao,8: 



<d.ud = S^±Re- 



±7 



(vud - ^')' ("'Pud'^s + ly 



(16) 

= ± Re,4^ ± 7 . f'^"-'... , (17) 



■(^2_^2)2 ' (a^2^<^^ + 1)2 



{aifl^ifs + 1)2 



(18) 
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